Invariant Perturbation Theory of Adiabatic Process 
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In this paper we present an invariant perturbation theory of adiabatic process according to the 
concepts of [/(l)-invariant adiabatic orbit and [7(l)-invariant adiabatic expansion. The probabilities 
of keeping the adiabatic orbit in the first-order and the second-order approximation are calculated, 
respectively. We also give a convenient sufficient condition. 
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Quantum adiabatic theorem and adiabatic approxi- 
mation are important contents in quantum mechanics 
P, [3, S Since the establishment of the quantum 
adiabatic theorem it has been broadly applied both in 
theory and experiments. In the deep investigation into 
quantum adiabatic process a lot of important results 
have been obtained, for example, Landau-Zener Tran- 
sition 0, Gell-Mann-Low theorem [6], Berry phase Q, 
holomonyQ. Recently, the content of quantum informa- 
tion and quantum computation has revived interest in 
quantum adiabatic theorem and adiabatic approximation 
[9, 10; 11, 12, 13, 14]. More recently, however, the con- 
sistency of the traditional adiabatic approximation con- 
dition has been doubted llSl 11611 and some new adiabatic 
conditions are proposed [l7l llSl. [2^ . which have caused 
confusions. Then, it is important to find a proper and 
convenient condition under which the evolution of the 
system can be considered adiabatic. 

This paper is aimed to propose a general invariant per- 
turbation theory for the quantum adiabatic process in 
terms of the concepts of [/(l)-invariant adiabatic orbit 
and J7(l)-invariant adiabatic expansion. Then, the prob- 
abilities of keeping in adiabatic orbit are obtained in the 
first-order and the second-order approximation, respec- 
tively. At last, a new and convenient adiabatic condition 
is presented. 

Firstly, we explain the invariant adiabatic basis and 
invariant adiabatic expansion with time-dependent 
coefficients. Let us consider a time-dependent Hamilto- 
nian H{t) with the initial state |to,0) at time t = Q, 
here m denotes the initial value of dimensionless quan- 
tum number set. We introduce a dimensionless time pa- 
rameter T — E„i{Q)t/h and a dimensionless Hamiltonian 
/i(t) — H{t)/ Em{0), Em{0) is the energy of the initial 
state. The time-dependent Schrodinger equation reads 



= /i(r)|$,„(r)), |$„,(r))|,=o = |m,0) 



' dr 



(1) 



here T is time-ordered operator. 

Definition 1: A state |<i>„ (r)), is a dynamic evolution 
orbit of system when the state |$„ (r)) satisfies Eq.(IT|), 
describing an evolution orbit varying with time in the 
Hilbert space. 



Furthermore, if we consider r as a fixed parameter, 
we can always solve the following quasi-stationary state 
equation 



/i(t) \ip„ (r)) = e„ (r) \Lp„ (r)) 



(2) 



then we will have the adiabatic solution or adiabatic 
orbit \ipn (r)) and its corresponding eigenvalue e„ (r) — 

Of course, although with the same initial state \m, 0), 
the dynamic evolution orbit |$m (t)) does not coincide 
with the adiabatic orbit \ipm (t)), or they are not even 
close to each other. Furthermore, because of the Hermi- 
tian of h (r) , all of these adiabatic orbits form a complete 
basis of the system. 

We denote 7„m (''■) = ii'finW \ ^mW) and the dot 
means the differentiation with respect to time. An adi- 
abatic orbit multiplied by an arbitrary time-dependent 
phase factor still describes the same adiabatic orbit. 

Definition 2: The following adiabatic orbits 



(r)) = exp |^*y" [em(A) - 7™™ (T)]dA| \(P„i{t)) 

(3) 

are invariant, up to a time-dependent phase factor, un- 
der C/(l) time-dependent transformation of the adiabatic 
orbit 



(4) 



with /m(0) = 0. We define this adiabatic orbit with 
special choice of the time-dependent phase factor as f/(l) 
invariant adiabatic basis which keeps the initial value of 
dimensionless quantum number set m invariant. 

Definition 3: If the \(pm (t)) in Eq.([3]) satisfies follow- 
ing conditions 



{ipn (t) I ipm (t)) = , ynj^m 



(5) 



thenj^J^,*" (r)) is also the dynamic evolution solution of 
Eq.(lT|. We call this adiabatic orbit adiabatic evolution 
orbif^ of system. 

Generally speaking, in an arbitrary evolution process, 
the dynamic evolution orbit \^m{T)) starting from the 
initial state |to,0) will change or even vibrate rapidly 
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among some adiabatic orbits. This case can be described 
by the probability P staying in the adiabatic obit 



(6) 



Next our task is to find when the dynamic orbit is 
sufficiently close to the adiabatic orbit if Eq.© is not 
satisfied. Then, we will give the correct adiabatic ap- 
proximation conditions. 

In the U(l) — invariant adiabatic basis the dynamic 
evolution orbit reads 

|$„(r)) = ^c„(r)|$f-(r)), |<i>™(r)} |.=o - |m,0), 

n 

(7) 

with initial conditions Cm(0) = l,c„(0) — 0,Vn ^ to. 
The time-dependent coefficients Cfc(T) are governed by 



(8) 



where the diagonal elements of matrix M(t) are zero and 
the non-diagonal elements of M(r) read 



mf'^irMi^mfr), w^k" 

OT 

d'n{ek'{ri)-ek''{ri))+^k'k"{ri), (9) 



here (r) = /([ d-q (7„„ (77) (77)) + arg7,m„ (r). 

= Cmn IS referred as geometric potential of this sys- 
tem. And for geometric potential one can obtain further 
detailed analysis and application in our other papers [l^l . 

Then we can get the expanding coefficients with initial 
conditions (7(0) (c„(0) = l,cfc(0) = 0,Vfc m) 



C{t) = I^Texp 
The element of Eq.lHO]) is 
Ckij) = ( Texp 



I I M{X)dX 




I I M{\)dX 




C(0). (10) 



(11) 



Apparently, Eq. (fTT|) shows that the dynamic evolution 
is just an adiabatic evolution if Eq.([5]) is satisfied. In 
addition, since M(t) is Hermitian, the probability of the 
evolution is conservative, that is 



\Crn{T)\^ +J2\cki^ 



= 1, 



(12) 



which shows the time-dependent system considered is not 
a dissipative system. 

Secondly, we try to get the probability Pm (t) of keep- 
ing in adiabatic orbit |0„j(t)). In the time-dependent 
dynamic evolution process, the probability of keeping in 
the adiabatic orbit |<I>5^'°(t)), i.e., keeping the dimen- 
sionless quantum numbers invariant, is 



Texp 



i I M(A) 




, (13) 



Then adiabatic approximation requires 



Prnir) 



(14) 



It means that the transition probability from dynamic 
evolution orbit to any other adiabatic orbits can be ne- 
glected. 

We can get the analytical expression of Pm{T) using 
the approach of coefficients ratio. Integrate Eq.® and 
get 



Cn(T) 



M(A)„; 



Cfc(T) 
Cm(r) 



Vn. 



(15) 



Obviously, this equation set can be solved by the iterative 
method. Consider Eq. (fTT|) we can get the probability of 
remaining in |$°,f"'(T)) from Eg. p^ 



Pmir)= n 

k^m 



exp 




jrexp 


i drjM (ri) 


^ km 


jrexp 


'iJ,^drjM{r^) 


^ mm 



M(\) 



mk 



(16) 



From Eg. pop we can also obtain a necessary and suffi- 



cient condition in a compact form 
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1 fe#m-^0 



dX 



(tcxp 






km 






) 


mm 



Mmk (A) 



0. 



(17) 



Thirdly, we give various approximate approaches to 
calculate the probability. 

1, According to Ea. ([T^ the first-order approximation of 

Prnir) is 



dx{<i>rix)\^TW) 



(18) 



Then the necessary and sufficient conditions can be de- 
scribed as 



dxin'^mnfw) 



0, Vfc m. 



(19) 



Eq. is of abundant contents which has been discussed 
in another paper 19]. Later we will obtain a convenient 
adiabatic condition from Eq. p^ . 

2, From Eq. (|lip the second-order approximation of 
PTnir) is 



Pmir) = 



1 - 



dXi 



dX2M (Xi) M {X2) 



(20) 

because the first-order term of Eg. pTj) is zero. The adi- 
abatic approximate reads 



E 

k^m 



dXi / dX2Mrak (Xl) Mk,n {X2] 

Jo 



0. (21) 



3, the first-order of the coefficients ratio is. Cfc(r), c„i(r) 
are all approximated in the first-order with Cm(r) = 1 



Cfc(r) 



Cm(T) jg 

We substitute Eq.jSH) into Eq.dTT]) 



dr/Mkmiv)- 



I k^7n ^ ^ J 



(22) 



(23) 



It should be pointed out that this method can only be 
applied to the situation without reversion of quantum 
state which can not be treated as adiabatic process with 
c™(t') = 

Fourthly, as application of adiabatic invariant pertur- 
bation theory we consider a spin-1/2 charged particle in 
magnetic field. Eq. (fT9|) can be rewritten as 
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—^0, Vfc 7^ m 
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Jo 





(24) 



Now we suppose q;(t) = 0, then a^r) can be linearly 
expanded 



(25) 



where ao, ^0 are constants. \^km (A)| is a periodic func- 
tion can be expanded in Fourier series 



dAe*"(^) l7fcm(A)| 



-,(fcm) 



(ilo + ^km, l) 



(26) 



The adiabatic approximate condition reads 



{flo + flkm.j) 



< 1, ci(T)=0, 



(27) 



which is a convenient sufficient conditions j2l| . In fact, 
Q!(r) = is a physical requirement. For this kind of 
Hamiltonian |22l | 

Hvit) - e-'*^He'*^, (28) 

where V and H — J2n -^n\En){En\ are two arbitrary 
time-independent Hamiltonian with {En, \En)} being the 
eigensystem of H. The quasi-stationary state of Hv{t) 
is e~**^|£'„). The adiabatic orbit reads 

l^m'M) = e-''^"^+''<-^"l^l^">^e-'*^|£;„). (29) 
Then we get the elements of the matrix M 



= ($r*(r)|$;^*(r)) 

^ g-«(E™-i5„)r+i(E„|V|iJ„>r-»(i5,.|V|£;„>r^^^jy|£,^ 



(30) 



The terms about t in the exponent are all linear, which 
means a = is satisfied. Therefore, we can prove 
concisely that the adiabatic condition Eq. ((27l) is suffi- 
cient for this kind of Hamiltonian. Next we will ex- 
plain our above adiabatic condition Eq. (|77|) with model 
in [16]. For system a, the adiabatic approximate con- 



ditions is 



^ 1, under this condition the 



LJQ+UJ COS t 

probability of remaining in adiabatic orbit is P„ = 
1 - ^iilpiisin^^ 1; For system 6, the adiabatic 
approximate conditions is |tan0| <C 1, under this con- 
dition the probability of remaining in adiabatic orbit is 
Pm = 1 — sin^ 9 s\T? ^ — ^ 1. It is not difficult to see that 



4 



the adiabatic conditions for system a and system h are 
different, then the inconsistency showed by this example 
does not exist when applying our condition. The trivial 
and special example of our condition with constant ge- 
ometric potential and energy gap is proved in [18| . Of 
course, our condition with time-dependent geometric po- 
tential is more natural and general. 

In conclusion, we present the invariant perturbation 
theory based on the concepts of invariant adiabatic orbit 
and expansion stated in our paper. The probability of 
keeping int the adiabatic orbit is given. Furthermore, 
we give an convenient adiabatic approximation condi- 
tion which is more convenient to apply than the general 
sufficient showed in [Im. The derivation of our condi- 



tion is more concise and general than [2^. The result 
is that we can obtain the approximate dynamic solution 
in perturbation theory by using the corresponding quasi- 
stationary equations and energy for nondissipative and 
smooth time-dependent process. 
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